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Abstract. Isotope shifts and hyperfine structure of tungsten were studied in the near 
UV range. We have used laser induced fluorescence spectroscopy on a pulsed supersonic 
beam to probe the 5 Do — > 5 -Fi transition at 384.9 nm, 7 S3 — > 7 P4 transition at 400.9 
nm, and 7 S3 — >• 7 P3 transition at 407.4 nm. Three new magnetic hyperfine constants 
are reported for 7 P 3 , 7 P 4 , and 5 Fi states. The isotope shifts of the 384.9 nm transition 
are presented for the first time, and the isotope shifts of 400.9 nm and 407.4 nm 
transition are measured with an order of magnitude higher precision compared to the 
previous measurements. As a result, the nuclear parameters A and A re ; are extracted 
from the isotope shifts with an improved precision. 
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1. Introduction 

High resolution spectroscopic study of tungsten(W) is of interest in various fields of 
physics [U [21 El IH |6]. As one of 5d shell atoms, its isotope shifts and hyperfine 
structure allow enhanced understanding of nuclear structure in the deformed region on 
the nuclear chart. Biittgenbach j2] analyzed the hyperfine structure of 4d- and 5d-shell 
atoms to study contact interaction terms. Aufmuth [4] extracted nuclear parameter A 
from tungsten optical isotope shifts. 

More recently there has been great interest in using tungsten as plasma facing 
material in controlled fusion [5]. Particularly, the International Thermonuclear 
Experimental Reactor (ITER) plans to use W tiles on the divertor plasma [6], along 
with the spectroscopic diagnostic system [7j to estimate the tungsten influx rate for 
plasma edge modeling. According to Skinner [8], the 7 S3 — > 7 P4 transition of neutral 
tungsten at 400.9 nm is considered promising due to its high transition probability. 
However, there is a complication in spectroscopic diagnostics due to singly ionized 
tungsten lines(W II) coincident ally at nearly the same wavelength as the neutral lines(W 
I), causing line blending issues [Bj. Another problem is the incomplete information on 
7 S3 — > 7 transition, as the latest study did not reveal the hyperfine structure due to 
limited resolution [4]. This would lead to further difficulties in accurate modeling of the 
tungsten influx rate at ITER. 

In this paper, we present isotope shifts and hyperfine structure measurements of 
three neutral tungsten transitions in the near UV range. Laser Induced Fluorescence 
(LIF) spectroscopy on an atomic beam were performed on 5 D — > 5 i*\ transition at 384.9 
nm, 7 S3 — > 7 transition at 400.9 nm and, 7 S3 — > 7 P3 transition at 407.4 nm. Magnetic 
dipole hyperfine constant of 7 S3 level is compared with the prior report [9], and three 
new measurements are made for the constants of 7 P3, 7 P4 and 5 F± levels. 

The isotope shifts of 384.9 nm transition is presented for the first time and the 
isotope shifts of 400.9 nm and 407.4 nm were measured with approximately 10 times 
higher precision compared with the previous results. High resolution isotope shift study 
of these three transitions enabled us to extract the nuclear parameters A and A re /, 
which were compared with the previous results [I]. Also, completely resolved tungsten 
spectrum at 400.9 nm is directly related to the application at ITER [8j. Experimental 
details are shown in section 2 and 3, followed by analysis in section 4,5 and conclusion 
in section 6. 

2. Experimental Methods 

The beam production of refractory metals has been limited due to its high melting 
points. Having the highest melting point of all metals, tungsten atomic beam cannot 
be generated efficiently through conventional oven method [TO]. In the past, we have 
used a resistively heated tungsten wire to generate a high flux beam [TT] . However, 
the Blackbody radiation coming from the heated wire was limiting our ability to detect 
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Figure 1. Diagram of Tungsten Beam Apparatus 

fluorescing light from the atomic transitions. Here, we use Smalley type pulse supersonic 
beam technique [12] to overcome these issues. Only a brief review of the technique is 
given here and the details can be found in ref. [T2] . 

Tungsten atoms are ablated from a rod (American Elements, 99.9% purity) by 
the third harmonic of the Nd:YAG pulse laser (Quantel), while the solenoid gas valve 
(Parker, general valve series 999) entrains the atoms with Argon buffer gas. The atoms 
get cooled down through buffer gas collisions resulting narrow Doppler linewidth. The 
turbo pump with 1500L/s of pumping speed maintained the operating pressures to be, 
5 x 10~ 6 Torr inside the vacuum chamber. Diagram of the experimental apparatus is 
shown in Fig. [TJ 

A tunable cw Ti-SAPH laser (Coherent MBR110) generated the IR light, which was 
frequency doubled with LBO crystal inside a cavity (Coherent MBD200) to produce the 
light in 380 nm - 410 nm range. This wide range of tunability enabled us to probe all 
the transitions presented in this paper using the same laser system. The probe laser was 
focused with an intensity of ~ ImW/ cm 2 at the intersection point where the laser beam 
is crossing the atomic beam perpendicularly. The laser induced fluorescing light was 
collected by a spherical lens into a water cooled Photo Multiplying Tube(Hammamatsu) 
connected to the photon counter. The photon counts were recorded simultaneously as 
the wavelength meter(High Finesse WSU series) measures the frequency of the probe 
laser with an accuracy level of few MHz. From the methods that are shown in the 
Appendix section, we have estimated the total uncertainty of frequency measurements 
of, 1.7MHz for 384.9 nm transition, 2.1MHz for 400.9 nm transition, and 2.9MHz for 
407.4 nm transition. 

3. Experimental Results 

Three of the measured W optical transitions are shown in figure [2] with their relevant 
energy levels. The 384.9 nm transition is from the 5 D ground electronic state with 
the 5d 4 6s 2 configuration to the 5 Fi excited state with 5d 3 6s 2 6p configuration [13] . Two 
other transitions of 400.9 nm and 407.4 nm shares the same low lying metastable state 
of 7 S3 with 5d 5 6s configuration, and excited states of 7 P 4 and 7 P 3 both have 5d 4 6s6p 
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Figure 2. Tungsten transitions and related energy levels. 



configuration [T3] . 

Tungsten has four even isotopes of 180 W, 182 W, 184 W, 186 W, with nuclear spin of 0, 
and one odd isotope of 183 W, with nuclear spin of 1/2 giving rise to hyperfine structure. 
As shown in Fig. [3]-^[5j 400.9 nm, 407.4 nm, and 384.9 nm transitions have three, four, 
and two allowed electric dipole hyperfine transitions, respectively. 

Previous studies of the 400.9 nm and 407.4 nm transition reported 20 ~ 
30MHz precision on isotope shift measurements with the hyperfine structures being 
unresolved pE]. In this paper, we have achieved 10 times higher precision on the isotope 
shifts measurements of both transitions and the hyperfine structures are clearly resolved 
as well. These are shown in fig. |3]and|4j The ratio of the peak amplitudes agreed within 
10% level of the natural abundance ratio combined with the hyperfine line strengths, 
which was calculated from the Clebsch-Gordan coefficients. The small disagreement 
comes from the fluctuation of atomic beam intensity. The linewidth of the measurements 
were mainly limited by Doppler broadening, with a Gaussian linewidth of 5MHz. We 
extract the line positions of each isotopes from individual Gaussian fits, except for 
the merged peaks of 183 W(6) and 183 W(c) in fig. HI which we fit to two Gaussians 
simultaneously. The fit gave an frequency uncertainty of 200kHz in line positions which 
is negligible compared with the measurement uncertainty described in the previous 
section. 

As seen on figure El the isotope shifts of 384.9 nm transition were partly unresolved 
due to limited resolution, and the 180 W with 0.12% natural abundance is believed to 
be masked underneath the hyperfine state (b) of 183 W. As some of the peaks weren't 
resolved, we had to take a different approach to analyze the isotope shifts. First step 
was to take an average of multiple scans of the transition, in order to eliminate the beam 
fluctuation effect and get an exact line profile. After this, we used a multi parameter 
Gaussian fit to extract the isotope and hyperfine shifts. Based on natural abundance 
ratio and the intensity rule for the hyperfine transitions, we identified two peaks on 
the wing side in figure [5] as the hyperfine states of 183 W. Also, we noticed the center 
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Figure 3. 7 S3 — > ' P± tungsten transition measured by laser induced fluorescence 
spectroscopy. Dashed lines show the positions of the smaller peaks. 

of gravity of 183 W lying closer to the right side of the unresolved center peak, which 
indicates negative isotope shifts for the 384.9 nm transition. With the above constraints 
for peak assignments, the least square fit was performed to match the experimental data. 
The least square fitting curve is shown in red dashed line with a Gaussian linewidth 
of 5MHz which overlaps well with the experimental data shown in black hollow circles. 
The blue solid line is the result of the simulation with a linewidth of 0.5MHz to show the 
line positions, and the yellow dotted line is the center of gravity position of 183 W. The 
180 W isotope was excluded from the fit, as the experimental data doesn't have enough 
resolution to constrain its position. We note that regardless of the peak assignments, 
the upper limits of 20 MHz/amu for the isotope shifts and 80 MHz for the hyperfine 
splitting of the 5 Fi state, can be deduced from figure El No previous study of isotope 
shift exists for 384.9 nm transition, however, Gluck [H] has reported isotope shifts 
for this transition with a frequency resolution of 30 MHz, which is consistent with our 
result within their uncertainty of measurement. 



4. W Hyperfine Structure Analysis 

The hyperfine structure of 183 W isotope is only caused by magnetic dipole interaction 
as it has a nuclear spin of 1/2 [15]. In this case, the hyperfine energy levels are given 
by 

F(F + 1)- J(J + !)-/(/ + 1) 



E r 



Ej + hA- 



(1) 
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Figure 4. ' S3 — > "' P3 tungsten transition was measured by laser induced fluorescence 
spectroscopy. Dashed lines show the positions of the smaller peaks. 
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Figure 5. 5 D — >• 5 Fi tungsten ground state transition was measured by laser induced 
fluorescence spectroscopy. Black hollow circles are the experimental data, red dashed 
line is the least square fit, blue solid line is the simulated line positions, and the yellow 
dotted line is the center of gravity position for 183 W. 
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Table 1. The magnetic dipole hyperfine constant A was measured and compared for 
7 S3, 7 P3, 7 Pi, and 5 Fi electronic states of Tungsten. The 1 a errors are shown inside 
the parentheses in units of the last decimal quoted. 



Electronic 
States 


A (MHz) 


This work Ref. [9] 


5d 5 6s 7 S3 


505.5(4) 505.592(12) 


5d 4 6s6p 7 P 3 


496.2(6) 


5d 4 6s6p 7 P 4 


440.5(6) 


hdHs 2 Qp 5 Fi 


47.1(11) 



where h is the Planck's constant and A is the magnetic hyperfine constant. The shifted 
energy levels due to this hyperfine interaction are shown on Fig. [3]-^ [5j The hyperfine 
frequency splitting within the electronic state becomes, 

^hyperfine = A (^J + . (2) 

It is straightforward to extract the magnetic hyperfine constant A from measured 
frequency splittings using equation [2J As some of the hyperfine splittings are measured 
in multiple transitions, such as the splitting of 7 S3 state, we take the weighted mean 
of the measurements from different transitions to get A, which are shown on Table [U 
We report three new hyperfine constants for 7 P 3 , 7 P 4 , 5 P 4 states, and the constant for 
7 S 3 state is compared with the previous result [9]. Good agreement is found within the 
error of the measurement for the 7 S3 state. 

Wyart suggested 5d 4 6s6p as the dominant configuration for 7 P 3 , 7 P 4 states, and 
5d 3 6s 2 6p configuration for the 5 Pi state [13]. Our measurements of larger constants A 
in the cases of 7 P 3 and 7 P 4 state compared with the case of 5 Pi state, support Wyart's 
configurations. This is due to the open shell 6s electron from 5d 4 6s6p configuration 
having bigger contributions to A than the 5d 3 6s 2 6p configuration. 



5. W Isotope Shift Analysis 

The theory of atomic isotope shifts are well developed in many publications [T6l [TTJ HU 
[19]. Due to the way our experiment is designed, we examine the total isotope shift of the 
electronic transition rather than the isotope shifts of individual electronic states. The 
total isotope shift of the electronic transition is equivalent to the difference in isotope 
shifts of upper and lower electronic states that are involved. Analyzing the isotope shift 
of individual electronic states will be discuss at the end of this section. 

Isotope shift between two isotopes with mass number A and A' of an optical 
transition i can be written as, 

° v i = du i,MS + d "i,FSi ( 3 ) 
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where MS stands for Mass Shift and FS is for Field Shift. The Mass shift is further 
separated into Normal Mass Shift (NMS) and Specific Mass Shift (SMS), both coming 
from the change in nuclear mass. 



° v iJiS = du i,NMS + du i,SMS ( 4 ) 

A' — A 

= {M^nms + M hS Ms) -jjr ( 5 ) 

The Field Shift (FS) originates from the change in volume and shape of the nucleus, 
therefore being directly related to the changes in mean square nuclear charge radii 5 (r 2 ). 
This relation is shown below, 

5v AA ' s = FiX AA> (6) 



X AA ' = 5 (rf A ' + g<5 (r*) AA ' + (r^' + ■■■ (7) 

^6(r 2 ) AA ' (8) 

Fi = Eif(Z), (9) 

where Ei, f{Z) are the electronic and relativistic correction factors defined in ref. [T9] 
and C n /C\ ratios are tabulated by Seltzer [T6] . 

The NMS can be calculated directly as shown by ref. [T7] . 

Mi N MS = ~ , (10) 
' t ,NMS lg361 ' V ) 

where z/j is the atomic transition in MHz. However, SMS needs to be extracted from the 
King plot [TS], which requires isotope shifts from two different transition, one of them 
being the reference transition with a known SMS. The ns 2 — nsnp transitions are often 
used as a reference transition because its SMS can be evaluated semi-empirically [17J. As 
we haven't measured any ns 2 — nsnp transitions in this work, we used the (5c? 4 6s 2 5 D\ - 
5<i 4 6s6p 1 Fi) 543.5 nm transition measured from W.G. Jin [2U] as a reference transition 
to make the King plots of our measured transitions. These King plots are shown in 
figure [6] where the axes are the modified isotope shifts of each transitions defined as, 

A A' 

^w = (^._^ MS )___. (11) 
From the linear fit, we obtain the relations, 

Slope = Ei , (12) 

-&543.5nm 

E- 

Intercept = M iyS MS ~ M U 3.5nm,SMS x ^ — - — • (13) 

-^543.5nm 
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Figure 6. King plot of the modified isotope shifts of 384.9 nm, 400.9 nm and 407.4 nm 
transitions with the 543.5 nm transition from W.G. Jin |20) as the reference transition. 
Factor of 184 ^ lg6 was multiplied to both axes for display purposes, which is to get 
the units in MHz. The 180 W isotope was not measured for 384.9 nm transition. 
Straight lines passing through the symbols are the least square fit results. Experimental 
uncertainties are all within the symbols. 



The SMS of 543.5 nm transition is estimated by a semi-empirical relation, 

^ U 543.5nm,SMS = (0 ± 0-5)SlS A ^ 5nmNMS , (14) 

which is only valid for ns 2 — nsnp transitions [17] . Combining equations 12 ~ 14, we 
can calculate the SMS of our measured transitions. Knowing NMS and SMS, the FS is 
readily calculated from eq. |3]and[5j The extracted values of NMS, SMS and the FS for 
384.5 nm, 400.9 nm and 407.4 nm transitions are given in Table |2j As discuss previously, 
isotope shifts of the 384.9 nm transition are given by the least square Gaussian fitting 
method. 

From eq. |6] and |9j the Field Shift is linked to the nuclear parameter X AA ' through 
Ei and f(Z). The electronic factors E 3 sA.9nm, E^Q^ nm and E i07Anm are derived from 
eq. [121 with an electronic factor E 5i3 5nm = 0.40882 calculated by Aufmuth |3|. The 
relativistic correction factor f(Z) was calculated by ref. [19] using the isotope shift 
constant of Blundell [21] . The extracted X AA ' from three different transitions are shown 
on the last column of Table |2j All the results show good agreement within the errors. 
Two independent sources of errors of X AA ' are shown inside the parentheses, the first 
parenthesis showing the experimental uncertainty from the field shift measurement and 
the second parenthesis showing theoretical uncertainty of Fi. The theoretical uncertainty 
in Fi comes from fractional error in the process of calculating E^ and f{Z), which is 
estimated to be 5% for all three transitions [22]. However, the experimental uncertainty 
of the field shift varies for different transitions. 

In order to see the individual error contributions during the process of FS derivation, 
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Table 2. NMS, SMS, FS, and X AA ' are shown for 384.9 nm, 400.9 nm and 407.4 
nm Tungsten transitions, ler errors are shown inside the parentheses. For \ AA , 
we show two independent errors, where the first parenthesis shows the experimental 
uncertainty from the field shift measurement and the second parenthesis showing 
theoretical uncertainty of Fi from equation |6] 



Ttp n ^it i on 

J. J. CXXlOl U 1VJ-L1 

(nm) 


luUliUUCO 

(A,A') 


(MHz) 


0U i,NMS 

(MHz) 


uu i,SMS 

(MHz) 


6v A £ 

uu i,FS 

(MHz) 


x aa> 
(fm 2 ) 


384.9 


184,186 


-12.7(1.6) 


24.8 


-60.2(0.9) 


22.6(1.9) 


0.0808(68)(57) 




182,184 


-10.0(1.6) 


25.3 


-61.5(0.9) 


26.1(1.9) 


0.0932(68)(65) 




182,183 


-4.5(1.6) 


12.7 


-30.9(0.9) 


13.7(1.9) 


0.0489(68) (34) 


400.8 


184,186 


423.6(2.1) 


23.8 


-71.4(3.7) 


471.2(4.3) 


0.0810(7)(57) 




182,184 


498.4(2.1) 


24.3 


-72.9(3.8) 


547.0(4.3) 


0.0929(7)(65) 




182,183 


263.2(2.1) 


12.2 


-36.7(2.1) 


287.7(2.9) 


0.0489(5)(34) 




180,182 


334.0(2.1) 


24.9 


-74.6(3.9) 


383.7(4.4) 


0.0651(7)(46) 


407.4 


184,186 


463.6(2.9) 


23.4 


-95.5(4.1) 


535.7(5.0) 


0.0803(7)(56) 




182,184 


541.4(2.9) 


23.9 


-97.6(4.2) 


615.1(5.1) 


0.0932(8)(65) 




182,183 


286.9(2.9) 


12.0 


-49.1(2.4) 


323.9(3.8) 


0.0490(6)(34) 




180,182 


355.6(2.9) 


24.5 


-99.8(4.3) 


430.9(5.2) 


0.0653(8)(46) 



Table 3. The sources of uncertainties in Siy^pg 86 from different transitions of 
tungsten. 

Source of Uncertainty Transition (nm) 

384.9 400.9 407.4 



Statistical Uncertainty of j/ 184 ' 186 (MHz) 


1.6 


1.6 


1.6 


Systematic Uncertainty of ja 8 ,186 (MHz) 


0.06 


1.2 


2.4 


Uncertainty of AA A X Intercept of King plot (MHz) 


0.9 


1.1 


1.4 


Uncertainty of v^^^ m SMS x Slope of King Plot (MHz) 


0.2 


3.6 


3.9 



we combine eq. El [51 [12] and [131 to get, 

r AA' r A A' r A A' AA A a I 

5v i FS = 5u t -Su^ MS - — - - x Intercept- Su^ 5nmjSMS x Slope.(15) 

We analyze different sources of uncertainties in eq. [15]for all three transitions, which 
are shown on Table [3] for the case of the isotope shift between 184 W and 186 W. Note 
that the second term of the right hand side of eq. [HI is a well defined value from eq. [101 
therefore has no error. Row 1 through 3 of Table [3] shows uncertainties that are linked 
to experimental uncertainty of our frequency measurements. However the uncertainty 
described on the last row, which contains the SMS of the reference line, only rely on 
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the semi-empirical relation of eq. [TH This becomes the dominating error contribution 
in the case of 400.9 nm and 407.4 nm transition. Therefore, improving the experimental 
errors alone, which only affects first three rows of Table [3j would not further decrease 
the uncertainty of the FS of 400.9 nm and 407.4 nm transitions. Better estimation of 
the SMS is required for further improvement. 

Since the extracted X AA ' in Table El showed good agreement among different 
transitions, we take the weighted mean of these to report the final values of X AA ' and 
compare with the previous results of ref . [19] on Table HI The errors are presented in the 
same manner as before, showing the experimental uncertainty inside the first parenthesis 
and the theoretical uncertainty inside the second parenthesis. We notice experimental 
uncertainties being much smaller than the theoretical uncertainties for all isotope pairs. 

From the above discussion, it is clear that the limiting factors for high precision 
X AA ' measurement mainly comes from the uncertainty in estimation of Specific Mass 
Shift through semi-empirical relation, and the theoretical uncertainties in the Field 
Shift. First we discuss about the possibility of SMS calculation. Instead of using the 
King plot approach, there has been several attempts of multi-configuration calculation 
of SMS [23l EU EE] , however, as pointed out by Aufmuth [19] , reliable prediction was 
never made. Following the approach of ref. [25J, the specific mass effect perturbs the 
energy of the electronic state ip between mass A and A' by an amount of, 



E Pi ■ 

i>3 



4), (16) 



where p» is the momentum of the i th electron. Thus, the frequency difference in upper 
state of ip u and lower state of ipi can be written as, 

AE SMS {^ U ) - AEsms^) . 

ou u-+l,SMS — 7 ■ \ L <) 

This would be a more general way of estimating the SMS as it does not rely on the 
semi-empirical relation shown on eq. [HJ however, the knowledge on tungsten electron 
wavefunctions are required. 

Similarly, we can write down an expression for the field effect perturbing the energy 
of ip between mass A and A' by, 

AE FsW = ^m f{ Z)X^', (18) 

where |-?/>(0) 2 | is the non-relativistic electron charge density at the nucleus and f(Z) 
is the same relativistic correction factor which was mentioned above. The frequency 
difference between the upper state of if) u and lower state of ipi becomes, 

r aa' &EfsM - AEpsj^) 

0V u -+l,FS = 1 • ( 19 ) 

As indicated from eq. [TSJ the FS analysis could be used as a study of |-0(O) 2 |. This 
term could be used for calculation of electric fields inside the atoms, which are important 
in electron Electric Dipole Moment experiments 
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Table 4. X AA and X A ^ compared with previous results from ref. [19]. la errors 
are shown inside the parentheses. For our results, the first parenthesis shows the 
experimental uncertainty from the field shift measurement and the second parenthesis 
showing theoretical uncertainty of Ei and f(Z) from equation |9] 



Isotopes 
(A,A') 


Weighted X AA 


' (/m 2 ) 


Weighted X AA ' 


This work 


Ref. Q]5] 


This work 


Ref. [H] 


184,186 
182,184 
182,183 
180,182 


0.0807(5)(33) 
0.0931(5)(38) 
0.0488(4)(20) 
0.0652(5)(26) 


0.084(7) 
0.097(8) 
0.051(5) 
0.068(8) 


1 

1.1537(11) 
0.6047(22) 
0.8079(14) 


1 

1.154(4) 
0.607(5) 
0.808(23) 



The tungsten electronic wavefunctions have not been studied in detail so far, which 
makes the SMS and FS calculations described in eq. [16] ~ eq. [19] unavailable. At the 
current state, the best way to reduce the fractional error of the nuclear parameter is by 
defining a relative X AA ' as, 

■^rel = ^184,186 ' ( 20 ) 

The fractional error of A^ e f is much smaller than X AA ' , which can be seen on the 
last two columns of Table [4] The results are in very good agreement with ref. [19] , with 
3 ~ 16 times better precision. 

6. Conclusion 

We performed Laser Induced Fluorescence spectroscopy on a tungsten atomic beam to 
study isotope shifts and hyperfine structure of three different optical transitions. Three 
new magnetic hyperfine constants were shown, which agreed with Wyart's configuration 
of 5d 4 6s6p for 7 P 3 , 7 P 4 states, and 5d 3 6s 2 6p for the 5 F 1 state. The isotope shifts of 384.9 
nm, 400.9 nm, and 407.4 nm transitions were analyzed to give nuclear parameters A and 
X re i. We have shown the dominating error contributions in the process of Field Shift 
and the A determination. The X re i, which is an alternative way of representing the 
nuclear parameter, had smaller fractional error than the absolute A. Both cases showed 
good agreement with the previous studies as well as improved precision. The 400.9 nm 
transition is related to the spectroscopic diagnostics of tungsten influx rate at ITER. 

7. Appendix : Wavelength Meter Calibration 

The uncertainty in relative frequency measurement could be divided into statistical 
uncertainty and systematic uncertainty. First we describe the way we characterize the 
statistical uncertainty. We have measured one of the hyperfine splitting of Tungsten 
transition multiple times and compare the results with the previously known number. 
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Table 5. Our measurements of Yb isotopes and hyperfine structure are shown as 
relative shift from 176 Yb and compared with ref. [27] . The 1 a errors are shown 
inside the parentheses. As for our results, we show two separate sources of errors, first 
parenthesis showing the statistical uncertainties and the second parenthesis showing 
the systematic uncertainties in frequency measurements. 

Shift from 176 Yb (MHz) 
Transition 

This work Ref. [27] 



173 Yb (7/ 2^. 5 /2) -1426.3(1.6)(3.6) -1431.872(60) 

171 Yb (1/2 ^ 1/2) -1170.6(1. 6)(3.6) -1177.231(60) 

176 Yb (1.6)(3.6) (60) 

174 Yb 955.2(1.6)(3.6) 954.832(60) 

172 Yb 1952.6(1.6)(3.6) 1954.852(60) 

170 Yb 3235.5(1.6)(3.6) 3241.177(60) 

173 Yb (5/2 ^ 5/2) 3258.8(1.6)(3.6) 3266.243(60) 

168 Yb 4613.4(1.6)(3.6) 4609.960(60) 

171 Yb (3/2 ^ 1/2) 4758.3(1.6)(3.6) 4759.440(60) 

173 Yb (5/2 ^ 3/2) 4762.4(1.6)(3.6) 4762.110(60) 



Figure [7J shows a histogram of relative frequency shift measurements between 183 W(a) 
and 183 W(6) of 7 S3 — > 7 transition, which essentially gives the hyperfine splitting of the 
7 S3 state as shown on fig. [3J The results were compared with the previous measurements 
from ref. [9]. The width of the histogram corresponds to the statistical uncertainty of 
the measurement, which was given as 1.6MHz from the fit. 

In order to test for a possible systematic uncertainty in relative frequency 
measurements from the wavelength meter, we measured the known Ytterbium(Yb) 
isotope shifts. Ytterbium was chosen due to a notable high precision isotope shift 
measurements from Natarajan group [27]. 

Having relatively low melting point, Yb atomic beam was produced by conventional 
method of resistively heated oven, rather than the technique we described in section 2. 
Multiple apertures along the beam line minimized the Doppler linewidth. A tunable 
frequency doubled diode laser system (Toptica DL pro) generated 555.6 nm light for the 
LIF spectroscopy on 1 Sq — » 3 Pi transition of Yb. Yb has 7 isotopes, 168 Yb ~ 176 Yb, 
with 2 of them, 171 Yb and 173 Yb, having hyperfine structure. Detailed analysis on Yb 
spectrum is omitted, as it goes outside the subject of this section. 

We compare the center line positions of all of our Yb isotopes and hyperfine 
structure measurements with ref. [27], which are shown in Table [5j Figure El shows 
the plot of Table [5] in black squares, as well as the red solid line showing the least 
square linear fit. Ideally, this line should have slope of 1, and intercepting at 0. 
Our fit gave a slope of 0.9990(6), and the intercept at 2.0(1.7)MHz with 1 a errors 
inside the parenthesis. The error of the slope, which is 6 x 10 -4 , represents the 
fractional uncertainty over the measured frequency range. Accordingly, we can assign 
systematic uncertainty of 3.6MHz over the range of 6GHz in the case of Yb isotope shift 
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Figure 7. The histogram of relative frequency shift between 183 W(a) and 183 W(&) of 
7 S3 — > 7 Pi transition is shown for 53 independent measurements. The Gauss fit for the 
histogram is shown in red solid line and the reference line position is shown in yellow 
dashed line. The uncertainty of the reference measurement is 42kHz which is within 
the thickness of the dashed line. 
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Figure 8. Our relative frequency measurements of Yb isotopes and hyperfine 
structure are shown in black squares. The uncertainty of both ours and the reference 
measurements are within the symbols. The red solid line shows the least square linear 
fit of the data points. 



measurements at 555.6 nm transition. By assuming that the wavelength meter would 
perform in a similar way at the wavelength range of 385 nm - 410 nm, we can apply 
these results for the estimation of systematic uncertainty of tungsten transitions being 
studied in this work. 

In conclusion, we report 1.6MHz of statistical uncertainty and 6 x 10 -4 of fractional 
systematic uncertainty for the frequency measurements from our wavelength meter. The 
fractional systematic uncertainty converts to 60kHz for 384.9 nm transition, 1.2MHz 
for 400.9 nm transition, and 2.4MHz for 407.4 nm transition. Adding statistical and 
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systematic uncertainties in quadrature, the total uncertainty of frequency measurements 
are, 1.6MHz for 384.9 nm transition, 2.1MHz for 400.9 nm transition, and 2.9MHz for 
407.4 nm transition. 
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